Precision cosmological measurements: independent evidence for dark energy 
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Using recent precision measurements of cosmological paramters, we re-examine whether these 
observations alone, independent of type la supernova surveys, are sufficient to imply the existence 
of dark energy. We find that best measurements of the age of the universe to, the Hubble parameter 
Ho and the matter fraction Q, m strongly favor an equation of state defined by (w < —1/3). This 
result is consistent with the existence of a repulsive, acceleration-causing component of energy if the 
universe is nearly flat. 
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The current era in cosmology seems to be the first in 
which local astrophysical measurements are consistent 
with the generally accepted large scale cosmology. To 
provide some historical context, consider the period from 
1980 to roughly 1995. Inflation offered us a large scale 
model for cosmology, requiring fitotal = 1, which could 
not find verification in measurements on smaller scales. 
Attempts to dynamically determine fitotal (e.g., 0, 13]) 
consistently returned results of O to tal ~ 0.25 ±0.10. This 
led to the notion Q that, under the fitotal = 1 prior, 
there must be a bias between the distribution of light 
(e.g. galaxies) and mass (e.g. the dark matter compo- 
nent). Not only did the Univese have to be dark mat- 
ter dominated, the distribution of that dark matter had 
to be signficantly different than the the distribution of 
light. At the time, this was the only way to reconcile the 
small scale measurements with the large scale (inflation) 
requirement. 

In this note we reinvestigate whether recent determina- 
tions of cosmological parameters are sufficient, by them- 
selves, to imply the existence of dark energy - specif- 
ically, a component of energy with equation of state 
w = p/p < — 1/3. In the mid-90's several authors [1, Q 
analyzed aggregate data based on globular cluster ages, 
clustering of galaxies, big bang nucleosynthesis, and the 
Hubble constant and concluded that something like a cos- 
mological constant might be necessary to produce a flat 
Universe. However, the conclusions were not definitive 
at the time due to the large uncertainty in the observa- 
tional parameters. Our purpose is to update these earlier 
investigations, accounting for improvements in precision. 
We will argue that observations of key parameters such 
as the age of the universe to, the Hubble parameter Hq 
and the matter fraction Q m have become definitive in 
support of dark energy. One might question the need 
for this analysis in the post-WMAP era, but it is impor- 
tant to understand whether increasingly precise measure- 
ments are consistent with the concordance cosmology ob- 
tained from best fits of WMAP data. Indeed, given the 
dramatic nature and consequences of dark energy, it is 
important to understand the observational evidence for 
it as broadly and robustly as possible. 

Despite the impressive results of the type la supernova 



collaborations Q , it is still possible that dust Q , evolu- 
tion effects Q or exotic particle physics [i| might alter 
the interpretation of the extracted redshift-distance re- 
lation. For example, the axion models in 9] account for 
the dimness of distant supernovae by conversion of pho- 
tons into axions in background galactic magnetic fields, 
rather than through accelerated expansion. Exotic par- 
ticle physics models which are less well motivated than 
axions, but perhaps no more counterintuitive than the 
existence of dark energy itself, might in principle explain 
the supernova data without requiring acceleration. How- 
ever, the demonstration that a dominant component of 
energy with w = p/ p < — 1/3 is strongly favored by the 
observed values of cosmological parameters provides a 
direct and robust argument for acceleration. 

We seek evidence for a component which has equation 
of state w = p/ p < —1/3. Recall the Einstein equation 
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The sign of the acceleration R is determined by the sign 
of X)i(Pi + ^Pi), where the sum runs over all contribu- 
tions to the energy momentum tensor. Strictly speaking, 
w < — 1/3 is the threshold for a component to cause ac- 
celeration when it is the only form of energy. If other 
forms of energy are non-negligible the overall sign of the 
right hand side of ([I]) might still be negative (i.e., the 
universe is decelerating, albeit more slowly than other- 
wise) even in the presence of energy with w < —1/3. 
Asymptotically, though, the component with the small- 
est positive or most negative value of w will eventually 
dominate all others. We recall that a cosmological con- 
stant has w = — 1, while a dynamical scalar model with 
non-zero vacuum energy typically has — 1 < w < 0. Val- 
ues of w less than —1 violate the null energy condition, 
and are generally associated with instabilities [1(3] ■ 

Analysis of the 3 year WMAP data [ll| favors a nega- 
tive pressure equation of state for models with constant 
w when constraints on the matter energy density are in- 
cluded (i.e., from the Sloan Digital Sky Survey or the 
2dF Galaxy Redshift Survey). In this note we conduct a 
simpler analysis in which the priors are transparent and 
easy to state. 
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We find that best measurements of the age of the 
universe £o , the Hubble parameter H$ and the matter 
fraction Q m are sufficient to require the existence, dur- 
ing some cosmologically significant epoch, of a repulsive, 
acceleration-causing (w < —1/3) component of energy, 
assuming the universe is nearly flat. A relation between 
these quantities is obtained using Einstein's equation for 
a Friedmann-Robertson- Walker universe. The analysis 
itself is not necessarily new, but it can now be applied 
for the very first time with stringent constraints due to 
recent precision measurements of the relevant cosmolog- 
ical parameters. 

The age of the universe is given by 
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where we have taken w constant in time and neglected the 
radiation component as it is numerically small. We also 
assume flatness, which implies fide = 1 — ^m, and allows 
us to define the integral as I(£l m ,w). The quantities to, 
Hq and Q m then determine w. 

In the more general case, where the dark energy com- 
ponent has time varying equation of state w(t), the sec- 
ond term in the denominator of the integral in ^ (the 
dark energy term) is more complicated, having the form 
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If (1 + 3w(x')) > for all x < x' < 1, the dark en- 
ergy term (j4|) is always decreasing with increasing x, and 
the denominator in ([3]) is larger for all x than it would 
be in the special case w = —1/3, where ((4]) is constant. 
Therefore, if the dark energy never exhibits a repulsive 
equation of state, so w(t) > —1/3 at all times, the inte- 
gral is bounded above: 



I(n m ,w> -1/3) < 7(fi m ,-l/3) 



Similarly, we deduce 



I(Q m ,w>w*) < I(Q m ,w*) . 
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In other words, in the most general case, unless the dark 
energy behaved repulsively during some earlier epoch, 
the integral /, and hence the product toHo, is bounded 
above by I{VL m , —1/3). Using measured values of to, Ho 
and fi m , it is therefore possible to deduce that a repulsive 
epoch must have occurred. (Note an epoch with repulsive 
energy does not necessarily imply overall acceleration, as 
discussed.) 

We now review the best measurements of to, Hq and 
f2 m . Systematically combining the results of distinct 



measurements using different techniques, each with dif- 
ferent statistical and systematic errors, is challenging. 
However, our discussion at least allows a reasonable guess 
at current global best values and uncertainties for these 
quantities. Examples of more sophisticated Bayesian 
analysis are given in 

to'. Our approach is made possible by relatively re- 
cent measurements of to with unprecedented accuracy. 
In the past, estimates of to have been made by either 
using model-dependent estimates for the ages of globu- 
lar clusters or through nuclear cosmochronometry. The 
former method has traditionally suffered from the un- 
known role of convection and its effects on the lifetimes 
of low mass/low metallicity stars. Krauss and Chaboyer 
[IH performed a thorough Monte Carlo analysis that in- 
cludes these uncertainties, to arrive at a firm lower limit 
of 11.2 Gyrs for t . However, t as large as 15 Gyrs is 
still allowable. Using Thorium cosmochronometry, Sne- 
den and Cowan [l4[ also find a lower limit of 11 Gyrs 
for to but acknowledge that lower limit could range up- 
wards by another 3-4 Gyrs. For the reasons cited, we 
do not use these methods or observations in construction 
our argument for the most probable value of to. 

Improvements in the precision of measuring to have 
utilized the white dwarf cooling curve and Hubble Space 
Telescope measurements of the halo globular cluster M4. 
Measurements by Hansen et al. (2002) II 511 r eport a value 
of 12.7 ± 0.7 Gyr. Hansen et al. (2004)[ll update this 
age to 12.1 ± 0.9 Gyr. The major source of systematic 
error in this analysis involves estimating the lag time be- 
tween the age of the Universe and the formation of glob- 
ular clusters. Numerical simulations of the Milky Way 
and its globular cluster system by Kravtsov and Gnedin 
(2005) L J_7J indicate that the peak formation of Globular 
Clusters occurs at z — 3-5. Using a mean formation red- 
shift of z — 4 implies that Globular Clusters formed at 
1.2 Gyr after the onset of the Big Bang. This then leads 
to a lower limit of to = 12.4 Gyr and a mean value of 



to 



13.3] 



Gyr. 



H : For decades, measurements of H were plagued 
by noise and biased samples. Today, however, there is 
good reason to believe that we have a relatively precise 
measure for this parameter as well. The Hubble Space 
Telescope Key Project for determining the Cepheid Zero 
Point and subsequent distance determinations to nearby 
galaxies using the Cepheid Period-luminosity relation- 
ship have returned a value of 72 ± 3 km/s/Mpc [l8] . The 
major source of systematic uncertainty in that measure- 
ment lies in the distance to the Large Magellanic Cloud 
(LMC), to which the zeropoint of the Cepheid Luminos- 
ity scale is anchored. Freedman and Madore [l!| quote 
a total systematic error of ±7 km/s/Mpc, but recent im- 
proved distance estimates for the LMC (e.g., Benedict 
et al. (2002) [H and Sebo et al. (2l|) have served to 
lower this systematic error down to ±4 km/s/Mpc (see 
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Ngeow and Kangur (2006) [22j). Moreover, confidence in 
the precision of Hq, as anchored by the LMC distance, is 
reinforced by recent measurements that are completely 
independent of the distance to the LMC. In the past, 
these kinds of measurements were also available but they 
had sufficiently large random error that precluded them 
from providing meaninful constraints on the value of Hq 
as determined from traditional distance scale ladder tech- 
niques. The new observations are: 

1) Using a sample of 38 X-ray clusters in combina- 
tion with the Sunyaev-Zeldovich effect, Bonamente et al. 
(2006) HH derive a value of H = 77.6 ± 5 km/s/Mpc. 
While there may be systematics associated with the non- 
spherical shape of clusters, their sample may is suffi- 
ciently large (and much larger than past samples) that 
this problem is removed by averaging. 

2) Wang et al. (2006) [23] have examined a sample 
of 109 SN of type la and have discovered important new 
corrections for metallicity and absorption (by dust) in 
determining SN la peak luminosity. This recalibration 
leads to Hq = 72 ± 6 km/s/Mpc. An independent treat- 
ment of SN la has been compiled by Riess et al. (2005) 
(25| which yields a value of H = 73 ± 4 km/s/Mpc with 
possible systematic error of ±5 km/s/Mpc. 

3) Koopmans et al. (2003) : 2j| perform a detailed anal- 
ysis of a gravitational lens system (from which a direct 
determination of the distance can be determined using a 
model mass distribution of the lens) to find H = 75±6.5 
km/s/Mpc. 

Averaging these 5 different results together formally 
leads to 74 ± 2.5 km/s/Mpc (error in the mean). Direct 
averaging is crude, but gives a characterization of the 
uncertainty. Averaging over systematic errors as well, 
we assume Hq = 74 ± 5 km/s/Mpc in further analysis. 
In contrast, one could use only method 1 and 3 above (as 
they completely circumvent the LMC distance problem) 
to obtain 76 ± 6 as the relevant range. 

f2 m : In contrast, £l m remains the most weakly con- 
strained cosmological observable. There are two reli- 
able methods of measurement: dynamical determinations 
based on infall to clusters of galaxies and/or the nature 
of large scale structure (e.g., Bothun et al. [13]) or by 
fitting the Hubble diagram to distant objects. In the 
first case, an unbiased and fairly large sample is needed 
for precision; in the second case, accurate distance mea- 
surements of intermediate redshift galaxies are required, 
and such measurements are ultimately based on the su- 
pernova luminosity scale. In principle, O m is highly con- 
strained by the multi-parameter maximum likelihood fit 
to the WMAP data; but this is an indirect determination 
of il m (as well as to) In the spirit of this analysis, we seek 
to use values of f2 m that have been directly determined. 

Note, though, that O m is now usually determined by 
assuming a flat Universe as a prior constraint. For in- 
stance, a recent accurate determination of fl rn results 



from analysis of the power spectrum of galaxy cluster- 
ing. Assuming a flat Universe, Sanchez et al. (2006) [28| 
find 0„j = 0.237 ± 0.02. In addition, Mohayee and Tully 
(2005) [29j revisit the peculiar velocities of galaxies in the 
Local Supercluster to derive O m = 0.22 ±0.02. Schindlcr 
(2002) 30] summarizes all techniques to determine fl m 
(including the more unreliable approaches such as the X- 
ray cluster luminosity function, weak gravitational lens- 
ing, or galaxy cluster evolution) . That summary yields a 
modal value of Q m = 0.3 (which is likely a realistic upper 
limit given the WMAP model) but also shows that most 
large scale structure studies yield values of Q m in the 
range 0.20 - 0.25 (which is consistent with the work done 
in the 1980s). Averaging together the Sanchez et al. and 
Mohayee and Tully studies produces a well constrained 
value of f2 m = 0.23 ± 0.02. For discussion below we take 
a conservatively large range for fi m , assuming 0.15 — 0.25 
to be a one standard deviation range about the central 
value. 

Results: In Fig. [IJ we plot I(il m , w) for Sl m = .15, .20 
and .25. f2 m = .15 corresponds to the curve with the 
largest values of toH . Taking t = 12.4 Gyr and H = 
69 km/s/Mpc, which are each one standard deviation 
below the favored (central) values in our assumed error 
model, we obtain to-ffo = .9, which corresponds to the 
grey horizontal line in the figure. The implications can 
be read directly from the figure. If w was always greater 
than —1/3, then some or all of our parameters must be 
well below their central values. 

From Fig. [IJ we see that taking to, Ho and fl m to 
each be one standard deviation below their central value 
(so, t = 12.4 Gyr, H = 69 km/s/Mpc and Q m = .15), 
an epoch with w < —.4 or so is required, which is just 
negative enough to imply acceleration (i? > 0). Taking 
to = 12.4 Gyr and Q m = .15, one would have to, e.g., 
push Ho below 67 km/s/Mpc to have w > —1/3, and be- 
low 50 km/s/Mpc to have w > (no negative pressure). 

We compute the likelihood of no epoch with w < w* 
(for given w*) as follows. First, we assume uncorrelated 
Gaussian errors in all three parameters: to = 13.3 ± 1 
Gyr, H = 74 ± 5 km/s/Mpc and fl m = .2 ± .05 (all one 
standard deviation). That is, we assume that the proba- 
bility distribution for the actual value each of parameter 
is normal, with maximum at the central value and stan- 
dard deviation given by the error estimate. We then com- 
pute, for a particular value of w* , the total probability 
that the parameters take on values for which inequality 
(J6j) is satisfied. In practice, this was done using Monte 
Carlo. 

The results are displayed in Fig. [2] (top curve) . Using 
this error model the probability of no epoch with w < 
— 1/3 is less than 4 percent. This is an overestimate of 
the likelihood, since the model allows values of, e.g., to 
which are much too low: to = 12.4 Gyr is more plausibly 
interpreted a strict minimum than minus one standard 



deviation from the central value. Modifying the error 
model so that values of to < 12.4 Gyr are not allowed 
reduces the likelihood of no epoch with w < —1/3 to 
about 1.3 percent. This is represented by the middle 
curve in Fig. O Adding a similar constraint that Q m > 
.15 leads to the lowest curve in the figure, and a likelihood 
of no epoch with w < —1/3 of about 0.8 percent. Fig. [3] 
is identical to Fig. [2] except that we have increased the 
one standard deviation error for Hq to ± 7 km/s/Mpc; 
the existence of dark energy is still strongly favored. 

We conclude that, unless systematic errors are signif- 
icantly larger than currently recognized, best measure- 
ments of the age of the universe to , the Hubble parameter 
Ho and the matter fraction Q m strongly favor the exis- 
tence of a repulsive dominant energy component, also 
known as dark energy. These observations are indepen- 
dent of type la supernova surveys: specifically, they are 
not sensitive to uncertainties 0, d, @] which affect the 
direct measurement of the distance-redshift relation at 
large z. 
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FIG. 1: Curves in the w — to Ho plane, each of which is an 
upper bound on toHo, for Q m = -25, .20, .15. The allowed 
region is between the top and bottom curves, and above the 
horizontal line toHo = 0.9. This requires w less than —1/3. 
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FIG. 2: Probability that w was always greater than w* for 
a range of wi* and various cuts on to and fl m . See text for 
details. 
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FIG. 3: Same as Fig. [2]except with larger Hubble uncertainty: 
H = 74 ± 7 km/s/Mpc. 
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